NASA-CR-1 9371 3 


/&/Z2./ 

if 


FINAL TECHNICAL REPORT FOR: 


ACCELERATION OF POSITRONS IN 
SUPERNOVA SHOCKS 


Grant Number: NAG 5-1131 
Period: April 15, 1989 - April 15, 1992 


un 


fH 

O 

t/t 

<N 

r- 1 

m 

rr\ 

| 


r -4 


o 

CO 

o 

c 

r-4 

2 


o 



o 



O' 



\ 



fO 



o 


an 


re 



OF 
i n 

\ 0) 


U_ 



2 



O co 

CO -P 


* 

o */> 


V- o 



< o 

fT 


ad X 

• C 


u j co 

u «- 



a,— 


LU <t 

< o 


O > 

u 


u a 

IT fC 


< 2 

«H O 



oc 

LLI 


* r 


a «p •*-> 
m 'D u u 
lO O C 
r- cl 2 
n z d ^ 
f> m it: a. 


I on f- r\j ^ 

o' 2 ^ a 
u 3 u o 

| CC r-* ^ 

< i— c * 
</>*-• x • > 

<r on u u ■- 

zcoiac 

vdh o 


Principal Investigator: 

Donald C. Ellison 
Associate Professor 

Department of Physics, North Carolina State University 
Raleigh, NC 27695-8202 

919-515-7227, Internet: ellison@pyvax.physics.ncsu.edu 


North Carolina State University 
Raleigh, North Carolina 27695-8202 



SUMMARY OF PROJECT NAG 5-1131 


During this project we investigated the acceleration of leptons (electrons and 
positrons) in collisionless shock waves. In particular, we were interested in how 
leptons are accelerated in the blast waves existing in the remnants of supernova 
explosions. Supernova remnants (SNRs) have long been considered as the most 
likely source of galactic cosmic rays but no definite connection between SNRs and 
the cosmic rays seen at earth can be made. The most convincing evidence that 
SNRs produce energetic particles is the fact that they are observed to be radio 
sources. The synchrotron radio emission is almost certainly produced by relativistic 
electrons (and/or positrons) and is direct evidence that the SNR blast wave produces 
energetic particles. Only by understanding lepton acceleration in shocks can the 
rich SNR data base be properly used to understand cosmic ray origins. 

The favored mechanism for accelerating particles in shocks is Fermi accelera- 
tion whereby particles scatter back and forth across the shock and receive repeated 
energy boosts as they collide with the upstream and downstream plasma. This 
mechanism has received a great deal of study, however, virtually all of the theo- 
retical work has been directed toward understanding ion acceleration. The study 
of lepton interactions with shocks is much more difficult to treat analytically (and 
computationally) and has not received as much work. Our project was directed at 
the neglected aspects of lepton acceleration. 

We performed calculations using a Monte Carlo particle computer simulation 
which overcomes many of the problems which have prevented lepton studies in 
the past. We can treat thermal particle injection, transrelativistic effects, and the 
feedback of the accelerated particles on the shock structure in a direct way not 
possible with other current techniques. Below is a list of some of the important 
discoveries made during this project: 

1. We showed that the efficiency of lepton acceleration depended critically on the 
lepton injection energy. If it is assumed that leptons have a diffusion coefficient 
similar to protons, thermal leptons will not be accelerated by the shock. Only 
when the lepton injection energy is many times the thermal energy, will signif- 
icant numbers of leptons be accelerated. This comes about because the shock 
structure is modified and smoothed by protons, not leptons, due to their high 
mass and the smooth shock prevents low energy leptons from receiving strong 
energy boosts upon crossing the shock. We were able to quantify this effect 
and it is published in Ellison and Reynolds (1991b). 

2. We showed that, even when injection effects are not important, that proton 
and lepton distribution functions produced by shocks are quite different in the 
critical energy range for producing the observed synchrotron emission. Until 
this result, it was universally assumed that ions and leptons would have si mil ar 
spectra. We also showed that transrelativistic effects produced proton spectra 
that were not in agreement with standard results from radio observations, but 
that the lepton spectra were, in fact, consistent with observations (Ellison and 
Reynolds 1991b; Reynolds and Ellison 1992). 

3. We performed simulations of relativistic shocks (shocks where the flow speed is 
a sizable fraction of the speed of fight) and discovered some interesting effects. 
We first demonstrated the power of the Monte Carlo technique by determining 
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the shock jump conditions in relativistic shocks (Ellison and Reynolds 1991a). 
We were able to perform calculations including both protons and leptons that 
can not be done analytically allowing us to calculate the amount of electron 
heating in relativistic shocks not undergoing Fermi acceleration. We then pro- 
ceeded to determine how relativistic shocks accelerate particles. We started 
with a test-particle calculation (Ellison, Jones, and Reynolds 1990) which al- 
lowed us to determine the acceleration time. From there we proceed to a full 
nonlinear calculation and found that nonlinear relativistic shocks treat protons 
and leptons even more differently than nonrelativistic shocks. The transrela- 
tivistic effects on the shock structure from the heavy ion component reduces 
the lepton acceleration to a tiny fraction of the ion acceleration. This effect 
is dramatic even if high energy leptons (many times thermal energy) are in- 
jected, and was totally unexpected. Our results have important consequences 
for astrophysics! environments expected to harbor relativistic flows such as 
extra-galactic radio sources and accretion onto compact objects. This work 
was published in Ellison (1991a, b). 

Below is a list of the published work associated with this project. 
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